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Physics of glassy systems 
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In this talk I present some of the recent theoretical results that have been obtained on glassy systems like spin 
glasses or structural glasses. The physical principles at the basis of the theory are explained in a simple language 
(without using replicas) and the results are compared with large scale numerical simulations. Finally we introduce 
the generalized fluctuation dissipation relation that can be directly tested in experiments with the present day 
technology. 
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1. Introduction 

In this talk I will take the point of view 
that glassiness (roughly speaking the appearance, 
when decreasing the temperature, of a very large 
equilibration time, much longer of those that 
can be observed on human scale) is related to 
metastability and to the presence of many equi- 
librium states. 

Replica theory |I|g] is the most powerful tool to 
deal with systems with many equilibrium states. 
Here I will present the main results of replica the- 
ory stressing the basic ideas, which have been re- 
cently recognized to be stochastic stability and 
separability. 

In this talk I will first elaborate on the rela- 
tions among glassiness and metastability. Later 
on I will introduce some models of glassy systems 
(spin glasses, tilings and structural glasses). I will 
then present the general theoretical interpreta- 
tion of these phenomena, which can become more 
quantitative using the two principles of stochas- 
tic stability and separability and I will report on 
some numerical simulations which support these 
ideas. In the next section I will highlight the con- 
nection among this approach and the off equilib- 
rium fluctuation dissipation relations. Finally I 
will present my conclusions. 



2. Glassiness and metastability 

It is not simple to give a precise definition of 
glassiness. In this talk I will take the point of 
view that glassiness corresponds to the presence 
of metastability in an open region of parameter 
space, a new and unusual phenomenon. 

Let us recall the usual case in which we have 
metastability. We consider a system that under- 
goes a first order phase transition when we change 
a parameter. The simplest example is the ferro- 
magnetic Ising mode: the control parameter is the 
magnetic field h. At low temperature the magne- 
tization m(h) is given by m{h) = TOsSign(/i) for 
small h (?Tis being the spontaneous magnetiza- 
tion) : the magnetization changes discontinuously 
at /i = in the low temperature phase where 
rris / 0. 

If we slowly change the magnetic field from pos- 
itive to negative h, we enter in a metastable re- 
gion where the magnetization is positive, and the 
magnetic field is negative. The system remains in 
this metastable state for a quite large time, given 
by T{h) oc exp{A/\h\"), where a — d ~ 1. When 
the observation time is of order of T{h) the sys- 
tem suddenly jumps into the stable state. This 
phenomenon is quite common: generally speak- 
ing we always enter into a metastable state when 
we cross a first order phase transition by changing 
some parameters. 

We can also define a linear response suscepti- 
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bility which is given by 



3. Glassy models 



PXLR= lim VH*)a(0)}^ 



(1) 



If we start with the the state where to > at 
/i = and we add a positive magnetic field h at 
time 0, the linear response susceptibility is equal 
to 

d 

XLR = lim ^T7-m(t, h), (2) 

m{t, h) being the magnetization at time t. The 
linear response susceptibility is not equal to the 
equilibrium susceptibility that in this case is infi- 
nite: 
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(3) 



h=0 



This is the usual stuff that is described in 
books. We claim that in glassy system there is 
an open region in parameter space where, if we 
change the parameters of the system (e.g. the 
magnetic field h) by an amount Ah, we have that 
Xlr 7^ Xeq ■ In the case of spin glasses we expect 
that for \h\ < hc{T) we stay in the glassy phase 
{hc{T) increases when we decrease the temper- 
ature and there is a temperature T^ such that 
hciTc) = 0). 

In this region 



Am(t) = XLR^h 
Am{t) = XeqAh 



for 1<< i << T{Ah), 
for T{Ah) « t, (4) 



where in some cases one finds numerically that 
t(A/i) is has a power like behaviour (e.g. 
T{Ah) oc |A/i|-4). 

It is convenient to define Xirr as Xeq = Xlr + 
Xirr- The glassy phase is thus characterized by a 
non zero value of Xirr ■ If we observe the system 
for a time less that t(A/i) the behaviour of the 
system at a given point of the parameter space 
depend on the previous story of the system and 
a strong hysteresis effects are present. 

The aim of the the theoretical study of glasses 
is, at least from my point of view, to get a the- 
oretical understanding of these effects and to ar- 
rive to a qualitative and quantitative control of 
these systems. Before presenting these efforts it 
is convenient to describe some models which have 
a glassy behavior. 



Generally speaking glassy systems can be di- 
vided into two categories: 

• Systems with quenched disorder, e.g. spin 
glasses. 

• Systems without quenched disorder, which 
are often translational invariant systems, 
e.g. tilings and glass forming liquids. 

This distinction in two categories is important: 
although the two categories behave in the same 
way in the mean field approximation, it is quite 
possible that in finite dimension they have some 
different features Q. We will see later that the 
definition of an ensemble will be slightly different 
in the two classes of models. 

3.1. Spin glasses 

The simplest model of spin glasses [|| is de- 
scribed by an Hamiltonian 

Hui<7) = - ^ U{i, k)a{i)aik) - ^ ha{t), (5) 

i,k i 

where C/(i, k) is a quenched Z2 gauge field on the 
lattice at infinite temperature (in plain words the 
variables U are defined on the links, take the val- 
ues ±1 and are uncorrelated quenched random 
variables), a{i) is a matter field which also be- 
longs to Z2. We are interested to compute the 
average value of the free energy (and of the cor- 
relation functions) at a temperature /3~^ of the 
matter field: 



F{i3) = -rMzum 

Z[;(/3)=^exp(-/3i/[;(a)). 



(6) 



{<y} 
A few comments are in order: 

• The overbar denote the average over the U: 
we are in the quenched approximation. The 
unquenched case in not physically interest- 
ing. 

• The quantity h is the magnetic field, which 
breaks gauge invariance. 



• At h — "to find the minimuni oi Hij(ay^ is 
equal "to find the Landau gauge", a being 
the gauge fixing. For given U in more than 
two dimension the problem of finding the 
global minimum of Hu{a) is an NP com- 
plete problem. 

• Gribov ambiguity (which follows from the 
NP completeness of the problem) implies 
that there are many local minima of Hu (a) 
(i.e. configurations whose energy does not 
decrease when we fiip one spin); their num- 
ber increase exponentially with the volume. 

As we have already remarked there is a a glassy 
region for not too large magnetic fields at low 
temperature. 

3.2. Wang Tilings 

The definition of Wang tilings is quite simple 
|5) . In each point of the lattice there are variables 
a{i) which can take M different values and are 
called tiles. The Hamiltonian is given by 



^(^) = E E E,ia{i),a{i + f,)). 



(7) 
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The model depend on the functions £^^(i7, r), 
which we assume can take only the values and 
1. For each M the number of different models 
is 2^^^ , however some of them are related by 
symmetries. Particular examples of Wang tilings 
are the Ising model {M = 2), the p-state Potts 
models (M = p), the Baxter vertex models. 

If there are configurations such that H{a) = 0, 
we say that the model admits a perfectly matched 
tiling. 

In two dimensions for sufficiently large Al weird 
things may happen: 

• For M > 16 there are functions E^{(j,t) 
such that for open boundary there are 
global minima of H with H — 0, but they 
cannot be periodic. In a finite volume 
system with periodic boundary conditions 
there are no perfectly matched tilings. 

• For M > 56 there are functions E^{a,T) 
such that the problem of deciding if a con- 
figuration of tiles on a finite volume can be 




Figure 1. The 16 Wang tiles which lead to aperi- 
odic tiling: E^, {a, r) is zero only if the two faces 
in contact of two nearby tiles are equal. 



imbedded in an infinite volume ground state 
with _ff = is not decidable. In other words 
there is no computer program which can tell 
us in a a priori bounded time if a given 
configuration can be imbedded in an infi- 
nite volume perfectly matched tiling. We 
can trivially construct a computer program 
that stops only if the configuration cannot 
be imbedded in an infinite volume ground 
state, however if L is the size of the configu- 
ration, the maximum time before stopping 
increases faster than any computable func- 
tion of L. 

Not so much is known on the thcrmodynamical 
properties of those models Q : most of the inves- 
tigations deal with the structure of the ground 
state configurations. 



3.3. Glasses 

Simple models of glasses are particles of differ- 
ent M species with Hamiltonian 

^= E E E K,fc(2;aW-Xb(fc)),(8) 

a, 6=1, A/ i=l,N(a) k=l,N(b) 

where N{a) = Nc{a), J2a '^i^) — ^ ^^'^ ^^^ quan- 
tities c are the concentrations of each different 
species. The model depends on the M{M + l)/2 
functions Vafi{x) and on the M concentrations. 

In order to have a glass transition it is crucial 
that the model does not crystallize, i.e. that the 
crystal ground state energy is bigger than that of 
lowest energy amorphous structure. 

A well studied model of glasses is realized for 
M = 2, c(l) = .8, c(2) = .2 with a potential that 
has a simple Lennard- Jones form [pj. 

4. The theoretical interpretation 

One way to interprete the presence of 
metastable states is to assume that for large (but 
finite) systems the phase space of equilibrium con- 
figurations can be decomposed into many finite 
volume states or valleys (lumps, as suggested Ta- 
lagrand |^). 

In other word we have a breaking of the ergod- 
icity: the Gibbs measure can be approximately 
written the sum of smaller disconnect pieces. We 
warn the reader that we have to navigate between 
Scilla and Cariddi 0. We cannot directly take 
the infinite volume limit because it is quite likely 
that the correlation functions do not have any in- 
finite volume limit (at least in a naive sense). On 
the other hand, if we work on a very large, but 
finite, volume system, the notion of equilibrium 
state is physically intuitive, but it need a more so- 
phisticated mathematic definition than the usual 
notion of infinite volume pure states (such a defi- 
nition can be given and it can be proved that the 
picture I am presenting rigorously holds, at least 
in some long range models S). 

Finite volume equilibrium states are character- 
ized by the following properties: 

• Connected correlation functions in each 
state go to zero al large distances (cum 
grano salis). 



• The time to go from one state to an other 
state is exponentially large with the volume. 

In glassy systems one finds that these states 
at low temperature have the following unusual 
features 

• Chaos: when we change the parameters in 
the Hamiltonian of any finite amount (typ- 
ically as soon as we make a change greater 
than F~^/^, V being the volume of the sys- 
tem)) we have state crossing: stable states 
become metastable. 

• No Gibbs rule: states coexist in an open 
region of parameter space. 

The crucial step is the decomposition of the 
Gibbs measure for a finite system into states: 



(■) 



(9) 



where the weights w satisfy the relation J^a ""^a = 
1. Of course the previous formula is only approx- 
imate for a finite systems (there is a fraction of 
phase space with a small, but finite probability 
that cannot be classified into states) and should 
become more and more exact when the volume 
goes to infinity. 

It is usually assumed that all the states have 
similar intensive properties, for example the inter- 
nal energy, the magnetization, the linear response 
susceptibility do not depend on the state: 

iV-i(M)„ =m + 7V-i/2j„, 



Xc.=N-\M'r^ 



Xlr 



(10) 



In this case we can also define the so 
called Edward Anderson parameter qEA = 
A^~^ ^,-(((T(i))Q,)^, which should not depend on 
the state a. 

However if we compute the full expression for 
the susceptibility, we find that the 0(N^'^/'^) vari- 
ations in the magnetization do matter: the ex- 
pression for Xirr turns out to be: 



Xirr = '^WaW^(5a ~ 5^f ■ 

a, 7 



(11) 



The physical origine of this extra term is quite 
clear: when we increase the magnetic field, the 
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Figure 2. The function Pj{q) for one given 16^ 
sample. 



Figure 3. The function Pj{q) for an other 16^ 
sample. 



states with higher magnetization become more 
likely than the states with lower magnetization 
and this effect contributes to the increase in the 
magnetization. However the time to jump to a 
state to an other state is very high (it is strictly 
infinite in the infinite volume limit if non linear 
effects are neglected) and this effect produces the 
separation of time scales relevant for xlr and Xeq ■ 
Having understood that the coexistence of 
many states in an open region of parameter space 
with the aforementioned properties is a crucial 
feature of glassy systems, we face the problem of 
being more quantitative. This quill be done in 
the next section. 



5. A quantitative approach 

In order to be specific, let us consider the case 
of spin glasses. A very important object is the 
overlap among two configurations (e.g. a and r): 



q 



(a,r)==7V-i^a(z)r(z) 



(12) 



In a similar way we can define the overlap qa,-f 
among two states a and 7 as the overlap among 
two generic configurations belonging to the two 



states. We also have: 

qo.,^^N'^Y.^a{i))^{a{i)). 



(13) 



In a similar way we can define a generalized 
overlap: given a local operator A{i) we have: 



q^ 



N- 



'^E 



{A{^))^{A{^)), 



(14) 



If we take A{i) — <j{i) we get q^ — q (i.e. the 
usual overlap): if we take A{i) = H{i) we get 
q^ = q^ (i.e. the energy overlap). 

Let us consider a given finite system character- 
ized by some parameters that we call J. In order 
to describe the structure of its states, we should 
give all the Wa and all the q^^- Equivalently we 
could introduce the probability distribution of the 
overlap given by 

Pjiq) = 
Z-2 J2 eM-mi^)+H{T))S{q{a,T)~q){l5) 

« ^WaWj6{q - qa,t)- 

a, 7 

The structure of the states depends on the sam- 
ple: we show in figs. (Q)-(0) two Pj{q) for two dif- 
ferent samples of three dimensional spin glasses of 
size 16^. 
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Figure 4. The average value of the overlap distri- 
bution (dashed line) and the energy overlap dis- 
tribution (full line). 



Figure 5. The quantity 912934 — 2/3(912) ~ 
1/3^. 



Is is clear that this description (i.e. the function 
Pj{q)) depends on the sample and numerically 
the dependence does not disappear by increasing 
the size. In this case the only thing that makes 
sense is to introduce an ensemble and to describe 
the statistical properties of that ensemble. At this 
end there are a few possibilities: 

• To average over the quenched disorder at 
fixed number of spins N . 

• To average over A^ in a window. 

• To average over small quenched added dis- 
order for a given sample. 

It is evident that the first possibility is empty is 
no quenched disorder is present. 

The statistical description would amount to as- 
sign the probability V{{w}, {9}) of finding a given 
configuration of weights w and overlaps 9. Such 
a functional is not so easily parametrized, so that 
we can introduce its moments, e.g. 



P{q)=Pj{q), P(9l,92)=Pj(9l)^j(92). (16) 

It is evident that V{{w},{q}) contains the 
same information of the infinite set of functions 
(P(9),P(9i,92),P(9i,92,93). 



In order to put order into this mess and to 
restrict the form of all possible 'P{{w},{q} two 
principle can be assumed: stochastic stability and 
separability. 

5.1. Stochastic stability 

Roughly speaking stochastic stability implies 
that the ensemble we consider is a generic ran- 
dom ensemble and it does not have any peculiar 
property |§-|ll|. More technically it implies that 
if we add a random perturbation, i.e. we consider 
the Hamiltonian 



H{e)ia) = Hia) + eHR{a), 



(17) 



where Hf({a) is a random perturbation, the ex- 
pectation value of everything is smooth in e. A 
typical example of random Hamiltonian we can 
add is HR{a) = Ar"^/^Ez,fc^»,fccr(i)o-(fc), where 
A^ is the total number of spin (as usual) and Ri^k 
are random quantities with zero average and vari- 
ance 1. It is crucial in the argument that infi- 
nite range Hamiltonians are allowed as pertur- 
bation. Stochastic stability implies very strong 
constraints on T'{{w}, {9}): the simplest one can 
be written under the form of the so called Guerra 
relations: 

P{qi,q2) = lPiqi)Siqi - 92) + ip(9l)P(92).(18) 



These relations are very well satisfied in short 
range model os spin glass. Indeed they imply that 



2 9 ^ / 9 

912934 = 3 [Tu 



2 1 

4 



(19) 



This relation whi ch is w ell satisfied numerically 
(see fig. (||, where (712934 is a quantity of order one 
for T < .9). As we shall see in the next section 
stochastic stability has important consequences 
on the dynamics. 

Although for the moment we cannot prove 
that spin glasses and other random systems are 
stochastically stable, stochastic stability it is 
rather likely to be correct for system which do not 
have any symmetry (or if symmetries are present, 
it should hold if restricted to those observables 
which are left invariant by the action of the sym- 
metry group). 

5.2. Separability 

Separability is a more subtle property pl] , p^ . 
It implies that there is only one kind of significant 
overlap. For example in spin glasses it amounts 
to say all overlaps (depending on the operator A) 
are given functions of the spin overlap. In other 
words for each local operator A there is a function 
/ (9) (which may depends on the magnetic field, 
temperature . . . ) such that q^^^ — f^{qa.-y))- 

In the general case an infinite number of quan- 
tities (i.e. q^ N MA) characterizes the mutual re- 
lations among a and 7. If separability is assumed 
this infinite number is reduced to one. 

There are some arguments |1^ that suggest 
that separability implies the ultrametricity prop- 
erty: 



qan) > nnn(gQ_^),(j/3,-y) V/?. 



(20) 



If ultrametricity holds, the set of states of a given 
system may be ordered on a tree (as it is usual in 
taxonomy) and the form of the probability func- 
tional V{{w}, {q}) is very similar to the one found 
in the mean field approximation. 

Separability and ultrametricity may be less 
compulsory properties than stochastic stability: 
this point is at present under investigations, how- 
ever the validity of ultrametricity seems to be sup- 
ported by numerical simulations for spin glasses, 
at least in 4 dimensions itlSi. 
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Figure 6. The correlation function C{t,tw) in 3-d 
spin glasses in the low temperature phase. 



6. Off-equilibrium dynamics 

A very interesting question is what happens 
when we cool the system starting from a random 
(high temperature) configuration at time zero? 
The question is non trivial if, cooling the system, 
we cross a phase transition. The behaviour of 
the system in these conditions tell us something 
about the order parameter in the low tempera- 
ture phase (in the following we will discuss only 
the behaviour of an infinite volume system). One 
of the most spectacular phenomenon is ageing. 

Let us suppose that we cool very fast the sys- 
tem from an high temperature configuration at 
time 0. We can define a two times correlation 
function 



(^ [t, tw) 



1 ^ 
— ^ i(CTi(t^)(Ti(i,„ +t)) = q{t^,t^+t) (21) 



TV 



In the limit tu, — > 00 at fixed t we recover 
the usual equilibrium correlation function (i.e. 
limt^^oo C{t, tyj) — C{t) ). However in the region 
where tw and t are both large, the dependance on 
tw does not disappear: for example in fig. g we 
show the dependance of the correlation function 
on t in four dimensional spin glasses: as soon as 
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Figure 7. Three different form of the function 
P{q) and the related function S{q) 



length can be defined [|l5[|). The increase of the 
the dynamical correlation length ^(t) corresponds 
to rearrangements of the regions which belong to 
different phases (some shrink and some others ex- 
pand) and it affects the two-times correlation. If 
^{t) increases as power of time, the previous ar- 
guments suggest the validity of simple ageing. Of 
course when ^(i) hits the size of the box ageing 
ends, however this phenomenon never happens for 
an infinity volume system if we use a local dynam- 
ics. 

In a similar way we can consider a time depen- 
dent Hamiltonian 



H = Hn 



9{t - tw)'^hiai 



(23) 



hi being a random magnetic field (which in spin 
glasses is gauge equivalent to a magnetic field 
with constant sign). 

After cooling the system at time zero, we switch 
on the random field at time t^,. In this way we can 
define a two times response function (the precise 
name is relaxation function): 



l3S{t,t^) 



1 ^ 

-Y 



i dai{tyj +t) , 



(24) 



t is of order t^^ the value of t^^, matters. If simple 
aging holds, we have that C{t,tu,) ~ C{t/t^^) in 
this region. 

We notice that limt^oo linit„_>oo ^(t, i^) is a 
non zero quantity, (i.e. it is equal to qEA) while 
limt^_»oo limt_»oo C{t, tw) is zero at zero external 
magnetic field. The non commutativity of the two 
limits is a very clear signal of existence of more 
than one equilibrium state. 

Aging can be intuitively interpreted in the fol- 
lowing way. At large times the system is ordered 
in one given phase in regions of size ^(t), where 
^{t) goes to infinity when t goes to infinity. These 
regions move with time and their coalescence pro- 
duces the increase in the correlations length. It 
is natural to suppose that 



C{t,t^)^f{at + t^)/atra)) 



(22) 



(a similar formula was suggested to be valid also 
in infinite range systems where no correlation 



A great simplification comes out if we close our 
eyes on the actual time dependence of the corre- 
lation function and of the relaxation function and 
we plot S versus C |jl^,^. We can expect that 
such a plot has a well defined limit when i^ goes 
to infinity. In this plot we must distinguish two 
regions: 

• The region C > qea is the equilibrium re- 
gion t << tw Here the fluctuation dissipa- 
tion theorem implies that 



dS 
dC 



(25) 



• The region C < qea is the aging region 
where the fluctuation dissipation theorem 
cannot be anymore applied. Here we can 
define the function 



^<'-) - -f • 



(26) 
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Figure 8. Relaxation function versus correlation 
in the Edwards- Anderson (EA) model in Z? = 3 
T = 0.7 ~ jTc and theoretical predictions from 



eq. 
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Figure 9. Relaxation function versus correlation 
in the Edwards- Anderson (EA) model in 13 = 4 
T = 0.7 ~ jTc and theoretical predictions from 
eq. M) 11. 



which characterizes the relations among 
fluctuations and response in the aging 
regime. 

It is rather surprising to note that a dynami- 
cal version of stochastic stability Q implies the 
relation 



X{C) 



dqP{q), 



(27) 



which was conjectured to be valid in infinite range 
models ^^. 

There are three main kinds of the dynamical re- 
sponse S{G), that correspond to different shapes 
of the static function P{q) as shown in fig. (ph. 

• Case A correspond to no replica symmetry 
breaking and it is the usual case when two 
phases are present (e.g. ferromagnetism) : 
P{q) is a delta function. 

• Case B correspond to one step replica sym- 
metry breaking and it is supposed to be re- 
alized in structural glasses: P{q) is a the 
sum of two delta functions. 

• Case C correspond to breaking the replica 
symmetry in an continuous way and is sup- 



posed to be realized in glasses and maybe 
in tiling models: P{q) is a the sum of two 
delta functions plus a smooth function in 
between. 

There are many numerical results which con- 
firm the validity of this analysis and classifica- 
tion. Some results are presented in the next two 
sections. 

6.1. Numerical results for spin glasses 

In spin glasses it is possible using the paral- 
lel tempering technique to fully thermalize sys- 
tems up to O(IO^). System of this size could 
not be thermalized in a reasonable time (eg. let 
than 10^^ Monte Carlo sweep) using the conven- 
tional Monte Carlo technique. After thermaliz- 
ing hundredths of these systems one can measure 
the average function P{q) and using the relation 
( [2^ ) one can predict the form of X{C) and con- 
sequently the dependence of R on C. 

These predictions can be tested by taking a 
large sample (e.g. 0(10®) spins) and by measur- 
ing the correlations and the relaxation function 
in the off-equilibrium aging regime. These results 
can be seen in fig. (|^) for Z) = 3 and fig. ( |6.l[ ) for 
D = A. The result are quite impressive if we con- 
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Figure 10. Relaxation versus correlation in a bi- 
nary mixture. 



Figure 11. Temperature dependence of the pa- 
rameter m in a binary mixture. 



sider that no free parameter is present and that 
the two curves are obtained with a very different 
procedure. 

Similar results have been obtained for a tiling 
model M. 

6.2. Numerical results for structural 
glasses 

In the case of structural glasses analytic compu- 
tations |l9| and phenomenological considerations 
poll support the proposal that 



X{C) =m for C< C*, 
X{C) ^1 for C < C*, 

where 

m « T/Tc 



(28) 



(29) 



(i.e m = 1 at T = Tc, m = at T == 0). 
The quantity m/3 is roughly independent from the 
temperature and has the meaning of the effective 
temperature in the space of the valleys. In other 
words replica symmetry should be broken at one 
step level and the symmetry breaking parameter 
m should be one at the critical temperature and 
vanish linearly with the temperature. 

This kind of behaviour is exactly the one that 
has been observed in some long range model of 



spin glasses. Numerical simulations support the 
correctness of these predictions. For example we 
see in fig. 6.2 the numerical evaluation of the 



parameter m in a glass forming binary mixture. 



7. Conclusions 

We have argued that glassiness, metastability 
and the existence of many equilibrium states are 
phenomena which are strongly linked. The the- 
ory for studying these phenomena has been devel- 
oped in the framework of replica theory, however 
we can abstract from the usual replica approach 
two principles, stochastic stability and separabil- 
ity, which are enough to fully characterize the the- 
ory. 

The most impressive phenomenon is the pres- 
ence of a new form of fluctuation dissipation rela- 
tion in the aging regime, where the characteristic 
function X{C) is linked to a quantity that can be 
defined in the static, i.e the function P{q). Direct 
measurements of these fluctuation dissipation re- 
lations in real (not numerical) experiments both 
for glasses and spin glasses are actually under the 
way. The outcome will be a crucial test for this 
theoretical approach and I hope that some results 
will be available in the next future. 
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